Dynamics and thermodynamics of linear quantum open systems 
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We analyze the behavior of a network of quantum oscillators coupled with a number of external 
environments. We show that the dynamics is such that the quantum state of the network always 
obeys a local master equation with a simple analytical solution. We use this to study the emer- 
gence of thermodynamical laws in the stationary regime, achieved for sufficiently long times if the 
environments are dissipative. We show that the validity of the second law implies the impossi- 
bility of building a quantum refrigerator without moving parts (therefore, a quantum absorption 
refrigerators requires non-linearity as an crucial ingredient, as recently proposed by Kosloff and 
others citeKosloffl,Kosloff2). We also show that the third law imposes strong constraints on the 
low frequency behavior of the environmental spectral densities. 

PACS numbers: 03.65.Yz, 03.67.Bg, 37.10.Rs 



Understanding the origin of classical physics from 
quantum laws is important not only for fundamental but 
also for practical reasons In fact, more and more 
experimentally accesible systems manifest quantum ef- 
fects while having a large number of degrees of freedom. 
In such systems, many of which are being developed as 
part of the quest for the quantum computer [sj , it is of 
fundamental importance to study the emergence of ther- 
modynamical laws to understand the dynamics of 
heat transport and decoherence. We present here some 
general results concerning not only the evolution of lin- 
ear open quantum networks but also the emergence of 
thermodynamics in such systems. This may be relevant 
not only for quantum information processing and simu- 
lation [6] but also be helpful in understanding natural 
processes such as the high efficiency of energy transfer in 
light harvesting complexes in biological systems Q. 

We investigate the behavior of a network of oscillators 
coupled with external reservoirs as shown in Figure 1. 
This is a generalization of the Quantum Brownian Mo- 
tion (QBM) model, which is the paradigm for an open 
quantum system [l], lil-flo| . We present four main original 
results. First we show that the dynamics is such that: 
(i) the quantum state always satisfies a simple local mas- 
ter equation and (ii) such equation always has a simple 
analytical solution. Then we use these results to study 
the emergence of classical thermodynamical laws in the 
long time limit. Thus, we derive such laws from first 
principles in a quantum model by means of a straightfor- 
ward calculation. Our derivation enables us to show two 
new thermodynamical results: (hi) We show that it is 
not possible to create a quantum absorption refrigerator 
using linear networks. Such refrigerators are interesting 
systems that, having no movable parts, induce the en- 
ergy to flow away from a cold reservoir by coupling it 
with a system that is itself coupled with several other 
reservoirs. In the quantum regime, absorption refrigera- 
tors were recently proposed by Kosloff and others [lll.ll2| 



using a specific non-linear model. Our work shows that 
non-linearity is, in fact, an essential ingredient of a quan- 
tum absorption refrigerator. Finally, we show that: (iv) 
the validity of the third law (unatainability of zero tem- 
perature) requires a simple constraint on the type of en- 
vironmental spectral density, which for low frequencies 
must behave as uj a with a > 0. 




FIG. 1: We analyze the most general network of interacting 
oscillators coupled with bosonic reservoirs characterized by 
arbitrary spectral densities and initial thermal states. 

We consider a network of N oscillators (in D dimen- 
sions) coupled with external environments. Below, we 
analyze in detail models where the total Hamiltonian is 
Ht = Hs + He + Hi n t , where the system Hamiltonian is 
H s = P T P/2 + X T VX/2. The symbol X= [ Xl ,...,x K ) 
denotes a column vector storing the coordinates (K = 
ND). Similarly, P is a vector collecting all momenta (we 
consider unit masses and superscript T denotes the trans- 
pose). The K x K matrix V defines the couplings within 
the system. The environmental Hamiltonian is He = 

Y.i h e,i with H E ,i = Y,k( n k /2mk+mkU>k% Z 2 -* and 
the interaction is H int — Y^i i k^fk x i^k ■ ^ n what fol- 
lows, we will study the evolution of the reduced density 
matrix of the system, obtained by tracing out the envi- 
ronment as p — Tte(pse)- 

Before discussing results for models based on the above 
Hamiltonian it is worth presenting two results that are 
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valid for a larger class of models. In fact, under the single 
assumption that the evolution preserves the Gaussian na- 
ture of the quantum states (together with the hermiticity 
and trace of p) we can show two simple results: (i) p(t) 
always satisfies the following time-local master equation: 

P = -i[Ha{t),p] -iTij(t)[xi,{pj,p}] 

- iFij(t)\pi,{xj,p}] - Dij(t)[xi, [xj,p\] 

- Dij (t) [pi , [ Pj ,p]}- Fij (t) [xi , [p 3 ,p]\. (1) 

Here H R = P T M R 1 (t)P/2 + X T V R {t)X/2 + f T (t)X + 
f T (t)P is a renormalized Hamiltonian that includes 
time dependent renormalized couplings through Vn(t), a 
renormalized mass matrix Mr{£) and renormalized forces 
f(t) and f(t). The non-unitary part of the master equa- 
tion includes the effect of relaxation (through T(t) and 
r(i)) and diffusion (through the matrices D(t), D(t) and 
F(t)). Below, we give formulae for these matrices for 
the generalized QBM defined by the Hamiltonian de- 
scribed above. Moreover, we can show the following: 
(ii) The state satisfying equation (1) can be explicitly 
written as a function of time. For this it is conve- 
nient to use the characteristic function, %(a,i), which 
is the expectation value of the displacement: x( a : = 
Tr(p(i) exp(— i(Pk p — Xk x ))). Here a = (k x ,k p ) is a 
2i4f-component vector defining momentum and position 
displacements. We can show that x(a,i) is: 

X(a, t) = x(*(*)«, 0) exp(--a T £ F (*)<*) exp(iJ T (t)a). 

(2) 

In this equation the time dependence is hidden in the 
2K x 2K- matrices $(t) and £f(£) and in the 2K -vector 
J(t). This results are valid for all linear models (with 
quadratic Hamiltonians). In such cases (1) and (2) show 
that the state evolves by a combination of a phase space 
flow and a Gaussian modulation (as the environment only 
induces renormalization, friction and diffusion). 

Here we sketch the derivation of these results that, to 
the best of our knowledge, have not been noticed so far. 
The first result is a generalization of the master equation 
obtained for QBM in Q . The second one is a generaliza- 
tion of a well known classical result by Rieder, Lebowitz 
and Lieb (see Q for references) who obtained the explicit 
form of the classical attractor state for general harmonic 
networks (for QBM the analytic solution was discussed 
in [T3| ) • Interestingly, both results are valid for arbitrary 
linear systems and can be obtained analyzing the 
properties of the evolution super-operator for p(t) , which 
is such that p(t) = J(p(0)). Thus, it is simple to show 
that any super-operator preserving the Gaussian nature 
of quantum states (together with their hermiticity 
and trace) can be written in the position represen- 
tation as: J(z, z', t; zq, z' , 0) = (z\ J(\z' )(zq\)\z') — 

e i(x T biZ+x T b 2 Zo+Xo b 3 z +Xob4Z ) e i(cJx+c%Xo) x 

e -x T aiix-x T ai2Xo-xJa 22X o/ 27r det(6 3 ). Here we de- 
fined if -vectors Z = (z + z')/2 and \ — z ~ z ' ■ The 



matrices b m and a mn together with the vectors q 
parametrize any quantum evolution. They depend on 
the microscopic model, as described below. The way 
to show the validity of the master equation (1) from 
the propagator is simple following the method discussed 
in [l|, [9(. We can compute the time derivative of the 
operator and notice that all the dependence on the initial 
coordinates zq and z' can be eliminated by rewriting 
these terms linear combination of others which 

are proportional to the final coordinates z^ and to 
the derivatives of the propagator with respect to them. 
This implies the validity of the master equation (I). 
Derivation of the solution (2) is even simpler as it can 
be done by integrating over initial coordinates using the 
explicit Gaussian propagator. 

Coefficients appearing in the equation (I) and in the 
solution (2) are determined by the time dependent func- 
tions appearing in the propagator. They can be obtained 
from microscopic models such as the one described by 
the above Hamiltonian. We did this using the Feynman- 
Vernon method [13( to compute the propagator. It is 
simple to show that in this case only two properties of the 
environment determine the propagator: the spectral den- 
sity and the initial temperatures. The spectral density 



is 7(c) = E^ W H where Z«( W ) = £*C^C$*( W - 
ujk) /2mkOJk (these K x K matrices are diagonal only if 
each environment couples with a single site). We assume 
that the initial state of each environment is thermal with 
temperature XJ . The spectral density determines the dis- 
sipation kernel defined as 7(7-) = J Q duil(ui) cos(uit)/uj. 
The temperature dependent noise kernel is z/(r) = 
du}J2 l I^(u})(l + 2ni(uj))cos(ujT), where n;(w) is the 
expectation value of the number operator in the l-th en- 
vironment (i.e., l + 2n;(o;) = coth(w/2/csXi)) (we denote 
v{lo) — J2i coth(w/2fcsTi)). For this generalized 

QBM model all coefficients in (1) and (2) depend on the 
K x K matrix G(t), which is the solution of: 



(3) 



G + V R G + 2 / dr>y(t - r)G(r) = 0. 



with initial conditions G(0) = and G(0) = I (the renor- 
malized potential is defined as Vr = V — 27(0)). The 
Laplace transform of G(t), denoted as G(s), is G(s) = 
(s 2 I + V R + 2s%s))-\ where j(s) = / °° duo I \u) s 2 /{uj 2 + 
s 2 ) is the transform of the dissipation kernel. In this 
model the master equation (1) is such that M(t) = I, 
f(t) = = f(t) (i.e., there is no mass renormaliza- 
tion nor renormalized forces) and T(t) = D(t) = (i.e. 
there is only normal dissipation and no momentum dif- 
fusion). The renormalization and friction matrices are 

given as V R (t) = (GG^G - G)(G - GG^G)- 1 and 

2T(t) = (GG _1 G - G)(G - GG^G)- 1 . In turn, the 
exact solution (2) is such that J(t) = and $(£) is 



*(*) 



G(i) G(t) 
G(t) G(t)_ 



(4) 
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For dissipative environments G(t) decays to zero for long 
times. This is the condition that must be satisfied for 
the asymptotic state to become independent of the ini- 
tial state (as seen from (2), this happens when 
vanishes). In this asymptotic limit the state is deter- 
mined by the matrix T,p(t), which is such that Sf(^) = 



(0,0) 
T F 

,(1,0) 



(*) 



r (0,l) 
T F 

,(1,1)/ 



(*) 



Each K x K submatrix can be writ- 



, (n,m) 

ten as a F 



F (*)J 

(*) = Jo* /J dti<ft2GW(ii)i/(ti-t 2 )G( m )(i a ). 
Thus, depends on G(t) and on the temperature- 

dependent noise kernel. This matrix is related with the 
diffusive matrices appearing in the master equation. In 



fact, we can write D = Sym(VR<j£ > ' 1 '' -I- 2r<jW' - ' + a 

The above results enable us to study the long 
time regime analyzing the emergence of ther- 



(i,ih 



a 



n,m) 



modynamical laws. In this 
Re $™ dujLo n+m i m - n G{iuj)v(io)G(-iuj). We consider 
the network to be divided into non-overlapping regions 
Si (I = 1,...,R), each one of which is coupled with an 
environment Ei (generalization to multiple environments 
coupled with each region is discussed below). We 
study the energy transfer between the system and the 
environments by computing the time derivative of the 
expectation value of the system's Hamiltonian which 
can be written as d{{Hn))/dt = J2iQi- This equation 
is equivalent to the first law of thermodynamics, where 
the heat current entering into the system through the 



projector onto the 1-th region of the network. In the 
stationary limit, the conservation law Qi = is 

satisfied. Using the formula for a < p' m ' we find that the 
heat current is: 



Qi = 



k 



uduQi k (uj) coth(cj/2fc B T fc 



(5) 



where 



Q lk (tu) = ImTr(P Sl V r G(iLo)I<> k \u)G(-iu)), (6) 
= -2tt\ Aj {u)G{iu)A k {u)\ 2 j ^ k. (7) 

To derive the last equation we used the explicit form of 
G(s) together with the fact that every spectral density 
I®(u)) can be written as I^(lu) = Ai{u)AJ(uj). The K- 
vector Ai(ui) has non-vanishing components only for the 
sites coupled with the environment Ei . In the derivation 
we also used the condition Yte(u! , y(iu))) — ttI(lu) (which 
is nothing but the simplest form of the fluctuation dis- 
sipation theorem Q). From these results we find that 
the heat transfer matrix Qih{oS) has the following prop- 
erties: i) Qikioj) = Qu{u); ii) QifcM < if j ^ k; hi) 
Yli Qik{u) — 0. Moreover, we also find that the heat flow 
can be related with the coefficients of the master equation 
as Qi -)■ Tr(P Sl (D - 2rc4 M) )). Using the above results 
we can derive the following thermodynamical laws: 



a) Equilibrium: Energy flows through the system if and 
only if there is a temperature gradient. This is due to the 
fact that the heat flow into the j-th site can be written 
as 



Qt = -2 



> 

ujdujQik(uj)(ni(u) - rife(w)), (8) 



As Qa-(w) < 0, the heat current Qi vanishes if and only 
if the Ti = Tk for all I, k. This is equivalent to the 0-th 
law of thermodynamics: When all the environments have 
the same temperature then energy does not flow through 
the system and equilibrium is reached. 

b) Heat flows from the hot to the cold reservoir. This is 
nothing but Clausius version of the second law of thermo- 
dynamics and can be shown as follows: If 7} is the highest 
temperature, then ni(uj) > nfc(w) for all k. Therefore, as 
Qifc(w) < we find that Qi > 0. This implies that the 
hottest reservoir always injects energy in the system. Of 
course, the conservation law implies that the other reser- 
voirs absorb this energy. Heat currents are 



Q 



3« = ^E / 

k=£l J ° 



1-th region is Qi = Tr (p Sl V r a { p 1] \ , where P Sl is the have that Q 
. V .J . imnnssible t< 



wdw\Aj (uj)G(iuj)Ak(uj)\ 2 (m (uj)-n k (uj)). 

(9) 

It is worth noticing that this same argument can be ap- 
plied to show that the coldest reservoir always absorbs 
energy from the system. In fact, if TJ < T k for all k we 
< 0. For this reason, we find that it is 



impossible to build an absorption refrigerator using an 
harmonic network such as the one we considered. Thus, 
such type of refrigerator, recently discussed by Kosloff 



and co-workers [ill, |12| , requires non-linearity as an es- 



sential ingredient. This no-go theorem for the linear ab- 
sorption refrigerator is general as it applies for any spec- 
tral density of the reservoirs. 

There are assumptions in the above derivation which 
can be easily relaxed. Thus, we assumed that each re- 
gion of the system couples with a single environment. 
We can generalize this as follows: We assume that the 
l-th region couples with Ni environments labeled by the 
index ai (each having temperature T ai ). For simplicity, 
we also assume that all the spectral densities /C> ai ) arc 
the same (this can also be relaxed), i.e. I^ l ' ai \ui) = 
Nil"'(u))/Ni. In this case we can show that the heat 
flow into the region-Z is given by equation ([9]), where 
the temperature dependent factor (ni(uj) — n k (ui)) must 
be replaced by a factor (hi(u>) — hk(w)), where hi(oj) — 
J2 a . n ai (uj)/Ni is the average number of excitations of 
all the environments coupled with the region Si. In the 
high temperature limit, the temperature dependent fac- 
tor in ([9]) becomes proportional to the difference between 
the average temperature of the reservoirs, but this fac- 
tor depends nonlincarly on the temperatures otherwise. 
Thus, the second law in this case states that the heat 
always flows into the region of the system that is cou- 
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pled with the reservoirs with the largest average num- 
ber of excitations (or the largest average temperature in 
the high temperature limit). Equivalently, heat always 
flows away from the region coupled with the reservoirs 
with the smallest average number of excitations (tem- 
perature). Of course, the environment with the lowest 
temperature may be coupled with a region which is it- 
self coupled with hotter environments in such a way that 
the average number of excitations is not the lowest one. 
Then, heat would flow away from the such environments. 
However, it is clear that by enlarging the definition of 
'the system' and by using the previous argument, one 
can prove that the environment with the lowest temper- 
ature would absorb heat. The other assumption in the 
above derivation is that different regions Si do not over- 
lap. However, this can also be generalized using the fact 
that the obtained results are valid for arbitrary networks. 
Thus, the case of overlapping regions can be analyzed as 
follows: Consider a network with regions S2 C S\. If Si ,2 
couple with different environments -E1.2, this situation is 
physically equivalent to that of a network where Si cou- 
ples with E\ while S2 couples strongly with another net- 
work S' 2 which is itself coupled with E%. In such case, all 
previous results directly apply. As a consequence of this, 
the no-go theorem for quantum absorption refrigerators 
applies to all linear networks. 

Another version of second law is obtained defining the 
entropy flow into the system as S — J2i Ql/T, where T) 
is any monotonically increasing function of Tj . Then S = 
~C duj ^k Qik(u)(l/f k - l/f,)(n*(w) -n,(w)) < 0. 

The above thermodynamical laws are valid for envi- 
ronments with arbitrary spectral densities. Below, we 
demonstrate that the validity of the third law imposes a 
constraint on I{uj). For simplicity, we consider environ- 
ments coupled with single sites, with spectral densities 
I^)(uS) = "fjUj Pj 9(lo) (9{uj) is a cutoff function vanishing 
when ui > A and jj a coupling constant). When pj — 1 
the environment is ohmic while it is super-ohmic (sub- 
ohmic) for p > 1 (p < 1). Using this, we can show that: 
c) The third law of thermodynamics emerges if Pj > 0.. 
Thus, the heat flow into the j-th site is: 

/>oo 

Qj = / u 1+p > +Pk i 3 i k e 2 {u)dLu\G{iu) ]k \ 2 



E( 



auj/k B Tj _ e ~au/k B T k \ 



(10) 



For low temperatures the integral is dominated by low 
frequencies. Taylor expanding the integrand and assum- 
ing that all Tj are close to the average T we find: 

Qj = 7j7fcT 1+w+Pfc (Tj - T k )(2 + Pj + p k )a kj . (11) 
where ay = 47r(fc B ) 2 +w+^ |G(0) jfc | 2 r(2 + Pj + Pk )((2 + 



Pj +Pk)- Thus, for two reservoirs with the same spec- 
tral density (i.e., pa = Pb) and a small temperature 
difference we find that the heat flow is Qa = 2(p + 
l)"/ 2 a A BT 2p+1 AT. This implies that the entropy flow 
into site- A is S A = Qa/T a oc T 2p . Therefore, the en- 
tropy flow vanishes as a power law when T — > if p > 0, 
which is nothing but Nernst's version of the third law. 



Thus, we presented new results concerning the evo- 
lution of open networks of quantum oscillators: (i) we 
deduced an exact master equation and (ii) we obtained 
its explicit solution. We used this to show how thermo- 
dynamical laws emerge from first principles for general 
(ohmic and non-ohmic) reservoirs. We proved a no-go 
theorem for linear absorption refrigerators and showed 
that the validity of the third law imposes a constraint on 
the low frequency behavior of the spectral densities. Such 
constraint is not the same as the one reported in 12j . 
Thus, according to our results the validity of the third 
law is consistent with sub-ohmic, ohmic and super-ohmic 
environments whereas the sub-ohmic ones are excluded 
according to the fl2| . 
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